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and Sergey P. Suetin 

This paper is dedicated to the 60th Birthday of Francisco (Paco) Marcelldn 

Abstract. In 1986 J. Nuttall published in Constructive Approximation the 
paper I28| . where with his usual insight he studied the behavior of the denom- 
inators ("generalized Jacobi polynomials") and the remainders of the Pade 
approximants to a special class of algebraic functions with 3 branch points. 
25 years later we try to look at this problem from a modern perspective. On 
one hand, the generalized Jacobi polynomials constitute an instance of the so- 
called Heine-Stieltjes polynomials, i.e. they are solutions of linear ODE with 
polynomial coefficients. On the other, they satisfy complex orthogonality re- 
lations, and thus are suitable for the Riemann-Hilbert asymptotic analysis. 
Along with the names mentioned in the title, this paper features also a spe- 
cial appearance by Riemann surfaces, quadratic differentials, compact sets of 
minimal capacity, special functions and other characters. 



1. Pade approximants to algebraic functions 

John Nuttall, whose name appears in the title along with such distinguished, 
actually illustrious, colleagues, initiated the study of convergence of Fade approx- 
imants for multivalued analytic functions on the plane. Obviously, he was not the 
first to consider this problem; the best known result in this sense is a theorem of 
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Markov (or Markoff) |21j . see also |25j . which assures the locally uniform conver- 
gence of diagonal Pade approximants to Markov functions: if 

where cr is a positive measure compactly supported on M with an infinite number of 
points of increase, then the diagonal Pade approximants [n/n\ff (see the definition 
in Section [2]) to a, which coincide with the approximants of the Chebyshev or J- 
continued fraction for this function, converge to a uniformly on compact subsets of 
the complement to the convex hull of the support of cr, and the convergence holds 
with a geometric rate. This theorem applies in particular to functions as 

- (z^-l)V^-z. 

In the same vein, Dumas 113) studied the case of the function of the form 



/(z) = ((z - ai)(z - a2)(z - a3)(z - 04)) ' -z -\ ^ z, 

with points Oj € C in general position and the branch of the square root selected 
in such a way that / is bounded at infinity. Dumas observed that the poles of the 
Pade approximants to / can be dense in C. 

However, Nuttall was the first to abandon the real line completely and start 
a convergence theory in a truly complex situation. From the Dumas' work it was 
clear that in a general situation we can no longer expect uniform convergenc^ The 
appropriate notion is the convergence in capacity 30 , 311 142| . that was developed 
independently by Gonchar and Nuttall. This is an analogue of convergence in mea- 
sure, where the Lebesgue or plane measure is replaced by the logarithmic capacity 
of the set. 

Still, the question about the domain of convergence (even in capacity) remained: 
if the approximated function / has a multi-valued analytic continuation to C ex- 
cept for a finite number of branch points, then the single- valued Pade approximants 
[n/n]/ cannot converge to / in this whole domain. They must "choose" the appro- 
priate region of convergence where / is single-valued too, and the boundary of this 
region should attract a sufficient number of poles of [n/n\f. 

In [261 1271 I29j Nuttall generalized Markov's theorem by considering a class 
of hyperelliptic functions of the form ri -I- r2h~^^'^, where /i is a polynomial of 
even degree and simple poles, and rj are holomorphic functions (it was extended 
later to meromorphic functions in the work of Stahl |39| and Suetin |44| ) . For 
these functions he found the domain where the convergence takes place: it is a 
complement to a system of arcs determined by the location of the branch points. 
Nuttall characterized this set as having a minimal logarithmic capacity among 
all other systems of cuts making the approximated function single-valued in their 
complement (see e.g. Figure [l]). In [26j Nuttall conjectured also that this result is 
valid for any analytic function on C with a finite number of branch points. 

The complete proof of this conjecture was given, even in a greater generality, by 
H. Stahl in a series of papers |32l I33L I34L 1371 140] , under the only assumption that 
the singularities of the function / form a polar set, i.e. a set of logarithmic capacity 
zero, see Theorem |2.1| below. Stahl also characterized the analytic arcs forming 



^We cannot expect uniform convergence even along subsequences, as it was shown in |4l [5l 

[20] . 
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Figure 1. Poles of ttti and 7r72 for /(z) = (z + 1.2)-3/7(2_ 0.7- 
1.75i)i/'^(z - 1 - 0.8i)2/^ Clearly visible is a "spurious" pole (left) 
or a distortion of the location of the poles (right). 



the boundary of this domain as trajectories of a rational quadratic differential with 
poles at the singularities of /. They are also a case of the so-called Boutroux curves, 
see e.g. [3j. 

The general results of Nuttall and Stahl (and also of Gonchar and Rakhmanov 
|15L I16j ) deal essentially with convergence in capacity and weak (equivalently, n- 
th root) asymptotics of the denominators and residues of the Pade approximants. 
However, strong or Szego-type asymptotics is extremely interesting, at least in 
order to clarify the behavior of the spurious poles. As it was mentioned, poles 
that appear within the domain of convergence in capacity and that receive the 
name (coined by G. Baker in the 1960s) of spurious, floating or wandering poles 
|4H I44j . can become the main obstacle for the uniform convergence. For some 
classes of elliptic and hyperelliptic functions |17l I44L 145] . the dynamics of the 
spurious poles is completely determined by the properties of the Riemann surface 
underlying the approximated function. Moreover, in the elliptic case |45| there is 
only one wandering pole, which greatly simplifies the description of the asymptotic 
behavior of the Pade approximants. 

The results in 45j were obtained when the approximated functions could be 
represented as Cauchy integrals supported on the critical trajectory of certain qua- 
dratic differential with four poles when this trajectory consists of two disconnected 
components; the analogue for the Chebotarev set was treated in [2]. 

In this paper we analyze the strong asymptotics of the Pade denominators and 
the residues of the Pade approximants for a canonical class of algebraic functions - a 
generalization of the Jacobi weight on [—1, 1], revisiting and extending the results of 
J. Nuttall in |28] . We look at this problems from two different perspectives, which 
give us two formally distinct answers. One of the main goals is to understand the 
relation between these two asymptotic expressions, in order to get new insight into 
the nature of this problem and the methods we have used. 

At the final stage of the preparation of this manuscript we learned about the 
closely related work [1] where the authors also apply one of the techniques used 
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by us to the asymptotic analysis of the Pade denominators in a shghtly more gen- 
eral situation. Our works, although close in spirit, have a number of substantial 
differences, and thus are rather complementary than overlapping. 



2. Statement of the problem 

Let aj, j — be distinct and in general, non-coUinear points on the 

complex plane C, A :— {ai, . . . ap}. Let 2l(C \ A) denote the set of functions 
/ holomorphic at infinity and such that / can be continued analytically (as a 
multivalued function) to the whole C \ A diagonal Pade approximant to / G 
2l(C \ ^) is a rational function 7r„ = [n/n]f = Pn/Qn of type {n,n), that is, such 
that both Pn,Qn G Pn (where P„ denotes the class of algebraic polynomials of 
degree < n), which has a maximal order of contact with / at infinity: 

(2.1) /(z) -7r„(z) =0(1/22"+!) z-^oo. 

This condition may be impossible to satisfy, but following Frobenius, we can obtain 
the coefficients of _P„ and Qn as a solution to the linear system 

(2.2) i?„(z) :=Q„(z)/(z)-P„(z) =0(l/z"+i) as z ^ oo, Q„ ^ 0. 



Equations (2.2) form an undetermined homogeneous linear system. Although the 



solution of (2.2) is not unique, the Pade approximant (rational function) 7r„ — 
Pn/Qn is- Hereafter, (PmQn) will always stand for the unique pair of relatively 
prime polynomials determining 7r„, and Qn is taken monic. The degree of Qn could 
be strictly < n. Nevertheless, in the beginning we assume that n is a normal index. 



so that Qn is of degree exactly n and (2.1) holds 



The tight connection of the analytic theory of Pade approximants with the 
(complex) orthogonal polynomials is given by the fact that the denominators Q„ 
satisfy an orthogonality relation 

(2.3) <j>t''Qn{t)f{t)dt = {), fc-0,l,...,n-l, 

where we integrate along a closed Jordan curve encircling A. 

Let ^/ be the family of compacts K containing A and such that / has a 
holomorphic (single-value) continuation to C \ if . From the fundamental work of 
Stahl |34U35L [36l I38j it follows that there exists F e .ft/ of minimal capacity, called 
the Stahl's compact. It is comprised of a finite number of analytic curves that are 
trajectories of a closed quadratic differential (see |22) or |43| ) and has a connected 
complement in C. In the particular case of p = 3 this is a star-shaped compact set 
known as the Chebotarev compact, see Section |4j 

Theorem 2.1 (Stahl). (i) There exists a polynomial 

p-2 

Viz) = Y[iz-v,) 
fe=i 

such that the complex Green function for F is 

(2.4) G{z,oc):^ rJ^dt, Aiz) := f[{z ~ a,) , z e C \ F. 
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In other words, T is given by the level curves o/ReG(-,cxi) — that join zeros of A 
orV. Alternatively, T is made of the closure of critical trajectories of the quadratic 
differential —{V/A){z)dz'^. 

(ii) The normalized zero counting measure for Qn converges weakly to the equi- 
librium measure Ar of V , given by 



(2.5) dXr{z)^-J^dz, 

TTi y A{z) 

with an appropriate choice of the branch of the square root, 
(ii) With an appropriate normalization of Rn, 

— log |i?„(2:)| — Rc G(z, oo), n — >• c», 
where denotes convergence in capacity in C\T. 

Remark 2.2. Depending on the function / (and the corresponding class M.f), 
Stahl's compact is not necessarily connected. For instance, already for p = 4 it 
can be a tree (for f{z) = A^/'^{z)) or a union of two analytic arcs (when f{z) = 
A^/^(z)). In the class Af Stahl's compact however is completely characterized by 
its S -property, namely 

(ReG'(z,oo)) = - — (ReG(z,c3o)), z e r°, 



where n± are the normal vectors to r°. The ^'-property and its generalizations 
play a crucial role in many branches of analysis and mathematical physics, see 
e.g. [ll[3lfT6l[22l[24l. 

In this paper we concentrate on a canonical example of a function from 2l(C \ 
A). Namely, let aj £ j — 1, . . . ,p, be such that ai + ■ ■ ■ + ap = 0. We 

additionally assume that no proper subset of a^ 's adds up to an integer, so that the 
corresponding Stahl's compact T € is a continuum. This is a sufficient condition 
for an underlying Riemann surface having the maximal genus (see below). However, 
for more precise hypotheses, see Assumption |3.3| The simplest non-trivial example 
of this situation is when p = 3, when T is star-shaped (Chebotarev continuum). As 
it follows from the works of Stahl, such "stars" along with analytic curves are the 
main building blocks for a generic F. 

Let 

(2.6) /(z) = [](z-a,)"^-=l + ^^, 

j=i k=i ^ 

where the expansion is convergent in the neighborhood D := {z £ C : \z\ > 
maxj- \aj\} of infinity. We will agree in denoting by /^/^ the branch of the square 
root in O such that f^/^{oo) = 1. 

As it was mentioned above, our main goal is to find the strong asymptotics 
of the Fade denominators (5„, as n — > oo. We derive this asymptotics using two 
complementary methods. The first one, developed in Section |3j is based on the 
differential equation satisfied by / and is a combination of the original ideas of 
Nuttall from |28j with some new developments in the asymptotic theory of gener- 
alized Heun differential equations. The second method is the non-linear steepest 
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descent analysis of Deift and Zhou (sec e.g. [11 ) based on the matrix Riemann- 
Hilbert problem |14) solvable in terms of and _R„. In Section [s] we apply it 
exclusively to the case of p = 3. 

These two methods provide formally different expressions for the leading term of 
asymptotics of the Pade denominators and Pade residues. We find their comparison 
in Section [6] very illuminating. 

The form of the asymptotics for Qn and i?„ was actually conjectured by Nuttall 
in |27| in terms of a function solving certain scalar boundary value problem. We 
show also that our results match the Nuttall's conjecture, see Section [6] 



3. Heine and Stieltjes, or asymptotic analysis based on the 
Liuoville- Green approximation 



The key observation is that function / in (2.6 1 is semiclassical: it satisfies the 
ODE 

It can be proved by standard methods (see e.g. |28p that as a consequence, the 
Pade denominators Q„, the Pade numerators P„, and the remainders i?„ satisfy 
the Laguerre equations: 

Theorem 3.1. For each normal index n there exist polynomials hn(x) = x^^"^ + 
• • • e Pp_2 and Dn{x) = x^p^^ H e V2p-4, such that 

(3.1) Ahnv" + {A'hn - Ahl^ - Bhn) v' - n{n + \)D„y = 
is solved by Qnf and P„. 

Remark 3.2. In the case of Qn, the ODE is of the form 

Ahny'^ + {A'hn - Ah'n + Bhn) y'n - n{n + l)Dnyn = 0. 
Let us use the notation Zk^n for the zeros of the polynomials hn from Theorem 

p-2 

(3.2) hn{z) ='^{z - Zk,n)- 

fe=l 

In order to simplify the situation and concentrate on the main ideas we impose the 
following assumptions on the zeros of hn : 

Assumption 3.3. There exists a constant M > such that for all sufficiently large 

n, 

(3.3) \zk,n\<M and \AVnh'n{zk,n)\ > fc=l,...,p-2. 

In other words, all zeros Zk^n of /i„ belong to the disk \z\ < M , and they stay away 
from the zeros of AVn and from each other 



The general case requires the spherical normalization for h„ and has to be treated separately. 
We avoid further discussion of this situation for the sake of simplicity. 
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Observe that the second part of this assumption is completely innocent: in a 
general case, we will have a reduction in genus and in a number of cycles in the 
basis of the underlying Riemann surface, see below. 

From Assumption |3.3| it follows that for all sufficiently large n the zeros of Dn 
lie in the disk \z\ < 2M, see e.g. |22j . so that that the set 

e„ := {z e C : Ah„Dn = 0} 

is uniformly bounded. 
Define in D, 



where the branch is chosen such that Hn{z) — logz + 0(1) as z — )■ oo. It can be 
extended as an analytic and multivalued function to the whole C. 



The polynomial solution P„ of (3.1) is known as a Heine- Stieltjes polynomial, 



while the corresponding coefficient Dn in (3.1 1 is called a Van Vleck polynomial, see 



e.g. |22l I23j . Theorem 2.1 from [23] gives a global description of the trajectories 
of the quadratic differential In particular, it is a quasi-closed differ- 

ential with one trajectory emanating from each zero of A and ending at infinity. 
Combining techniques from ^23 and [28 we get 

Theorem 3.4. For any a £ A there exists a progressive patl^j — 7(a), starting at 
a point zq G D and returning back to zq, which is homotopic inC\ e„ to a contour 
7 with a G Int(7) and e„ \ {a} C Ext(7). 

For any such a progressive path 7 we have for z e 7 U D, 

(3.4) i?„(,)=a„,i^^(z)e-("+WW (l + 5i(z)). 

If for p > 0, dist(z,e„) > p, then n\5i{z)\ is uniformly bounded by a constant 
depending on p. 

Remark 3.5. This formula should be understood in the following way: the right 
hand side is chosen for z e D according to the branch of i?„ described above, and 
then both the left and the right hand sides are continued analytically along 7. In 
this way this formula may be extended from progressive paths to rectangles in the 
= e~^"'^^^ plane. 



Constant C = C„.i in (3.4) depends on the normalization of i?n 



Formula (3.4) is not totally satisfactory, since it has a number of undetermined 
parameters. Our next task is to clarify their behavior. 

Using Assumption [3]3] and compactness argument we can choose a subsequence 



A = {uk} C N such that 

hn — > h, Dn D, as n G A, 



so that by (3.4), 



n 




^ \og\Rn{z)\^-Re r J^^dt, neA, 



for z E O. Let F be the Stahl's compact associated with /, i.e. / is holomorphic in 
C \ F, and F has the minimal capacity in the class 2l(C \ A). Then, it follows from 

■^Progressive path 7 means that Rei?„(z) is non increasing along 7. 
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Theorem 2.1 that D/{Ah) = V/A, with V definmg F. Since A was an arbitrary 
convergent subsequence, we conclude that actuaUy 



hm<=hm-^ = ^ = G'(.,oo). 



This estabhshes 



Lemma 3.6. For the polynomials I?„ in (3.1 1 we have the representation Dn{z) 
Vn{z)hn{z), such that 

p-2 



Vn{z) ^ Y\{z - Vk,n) ^ V{z), rn>oo, 



fc=l 



and h.n{z) — z^ ^ + . . . satisfies ~ ^ as n oo. 



Observe also that from our Conjecture |3.3| it follows that all zeros of V are 
simple. 

Since iJ„ appears multiplied by n in (3.4 1, we need to estimate the rate of 
convergence of to G'. Together with the Green function G it is convenient to 
consider also 

G„(z, oo) := f\l^dt. 



Ait) 

The following lemma is an elementary observation: 
Lemma 3.7. We have 

H'Jz) = G'^{z, oo) (1 + 5h,n{z) + - G'{z, oo) (1 + ^,,,„(z) + SvA^) + ^2(2)) : 

where 

(3.5) Oh,n{z) ■■= ^, , s = > ^ , with Pk.n ■ 



(3.6) 5y,„(2;):=V' with Avk.n ■= Vk.n - Vk, 

^ — ' Z — Vu 

fe=l " 

e,^0{5lAz)) + 0{5lAz)), n^oo. 
Proof. Applying the identity ^Jl + i = 1 + + 0{(^), £, ^ 00, we get 



By Assumption |3.3[ zeros of /i„ are all simple for n large enough. Hence, using the 
partial fraction decomposition for (ft,„ — h„)/hn we obtain the second identity in 

Finally, differentiating GJj(-,oo) with respect to its parameters Vk,m we obtain 
that G;,(-,oo) = G'(-,oo) {l + Sv,n{^) + 0{dl.Jz))), with Sv,n given in □ 

In order to find the asymptotics for ~ 5h^n we need the following result: 

Lemma 3.8. At any zero z^^n of hn we have 

(3.7) n{n + l)Dl = Ah',D'^ - [Ah':, + Bh'^)D,, . 



A JOHN NUTTALL'S WORK 25 YEARS LATER 



9 



Proof. Differentiating (3.11 and evaluating the result at 2 = Zk,n we get 
(3.8) - {Ah'; + Bh'„ + N^Dn) K - N^KRn = 0, 

where 



(3.9) 



Also from (3.1), for z — Zk^n, 
(3.10) 



(3.8)-(3.10) give us a homogeneous linear system on (i?„(z), i?^(z)) with a non- 



trivial solution, since by the uniqueness theorem, at a regular point of ( |3.1[ ), both 
Rn and i?^ cannot vanish simultaneously. Hence, the determinant of this system 
is zero, which yields the assertion. □ 



Using that Dn = Vnhn we obtain from (3.7), 

N^V^hl = AKXVn + hn{AV:,h'^ - Ahl + Bh'^) for z = Zk, 
and since by Assumption 



3.3 





K - K 



Ah'l + Bh' 



As a consequence, we get the following lemma: 



for z = Zfe „. 



Lemma 3.9. For f3k,n defined in (3.5), 

1 ^Kn) (l+0(<5.,„))- 



1 A{zk^n) 



{l+0{8h,^)+0{5v,n))^ 



Observe that the last identity is obtained applying also Lem ma |3.7[ 
Next, we use the possibility of the analytic continuation in (3.4) in order to 
derive the asymptotic identities on the unknown parameters. 

Lemma 3.10. Let ^ be a cycle (simple closed curve) in C\e„ enclosing two points, 
say ai, 02 G A, in such a way that the rest of points from e^ are exterior to 7. Then 



(3.12) 
with 



N <b K{t)dt^T{^,f){l + 0{l/n)), 
'7 



T(7,/) :=±log 



SmTTCll 



+ 2nim, to G Z, 



sm 7ra2 

where the sign is uniquely determined by the branch of the square root and the 
orientation of the contour 7 chosen. 



Proof. By Theorem |3.4[ for any Oj € A there exists a progressive path 
from D to D that is a closed Jordan curve separating Oj from other points of e„; 
assume 7^ positively oriented with respect to aj . 

Observe that both analytic germs / and i?„ = Qnf — Pn in D allow for the 



■>7j 



Qnf' 



1] 



analytic continuations along any such a path. Denote by f^. and R 
the values of these functions that we obtain after the analytic continuations of / 
and i?„, respectively, along 7^. If we denote by — 7^ the negatively oriented contour 
7j , then for z S 7j H D , 
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Consider for instance paths 71 and 
Qnf - Pn, we have 



(3.13) 

But 

(3.14) 



"72 • By the definition of the residue, Rn 

f-r2~f, 



iz). 



f-r, - f 
fy.-f 



-27ria2 



1 



sm 7ra2 



7- 



71 



sm 7rQ;2 
sin TTOfi 



On the other hand, during the analytic continuation the residue of the Pade approx- 
imant picks up a dominant term; thus, i?„ is geometricaUy small in D in comparison 
with Rn,i, i = 1,2 (see (3.4|), and from Theorem 3.4 we have for z E O, 



Rn 



Rn 



Rn. 



71 



Rn 



(^) 



R„ 



R 



■n.l 



(3.15) 



^(z) expf(n + l/2) / 

/ 71 V 71 - 



-72 



Kit)dt] (i + o(i)) 



/7 



■-{z) exp TV 



i/;(i)dO (i + o(i)) 



(observe that the negative sign comes from the fact that the term A^^^'^ is multiplied 
by ±1 after its analytic continuation; orientations of 71 and —72 are opposite, so 
after division we gain the —1 factor). Identities (3.13|-(3.15) yield the assertion 
with 7 = 71 — 72 or any any cycle homotopic to it in C \ e„. In order to extend the 
theorem to an arbitrary cycle 7 in C \ e„ , we observe that if during the homotopic 
deformation of the contour we cross a pair of adjacent zeros of /i„ and /i„ (see 
Lemma 3.6), both Rn^-yi and R 
valid. 



gain a change of sign, so that (3.12) remains 

□ 



We introduce the Riemann surface TZ defined by the equation w'^ — A{z)V{z). 
It is a hyperelliptic Riemann surface that can be considered as a two-sheeted cov- 
ering oi C, TZ = {z = {z,w) € C^}, with two sheets, 7?,'^-' and TZ^'^\ cut along 
Stahl's compact F and glued together in the standard way. From Assumption 3.3 
and using the Riemann-Hurwitz formula we easily see that the genus of 7?, is p — 2. 
The canonical projection tt : 7?. — > C is given by tt{z) = z ioi z — {z,w) € C^. We 
denote z'^^^ — t:^^{z) n T?.^-'^ j ~ 1,2, and we convene that sheet TZ^^^^ over C \ F 
is specified by the condition w/z^ — > 1 as 2; — oo^-^^ € 7?.'^"'^^. In this way, function 
w = [AVy^'^ is single- valued on TZ. Unless specified otherwise, we identify the first 
sheet TZ^^'^ with the domain C \ F = tt{TZ^^''). We construct analogously the Rie- 
mann surface TZn defined by the equation w'^ — A{z)Vn{z) . Again, by Assumption 
|3.3[ the genus of TZn is p — 2 for n large enough. 

Note that a homology basis of cycles of TZ and TZn can be constructed from an 
integer combination of cycles jij = ^{ai,aj) = ji — "fj considered in the proof of 
Lemma 3.10 Thus, (3.12) is valid for any cycle 7 on 7?, with the right hand side 



r(7, /) depending on the representation of 7 in terms of the basis of cycles "fij. We 
select a homology basis of cycles Fj oi TZ, j — 1,2, . . . ,2p ~ 4, in such a way that, 
in the standard terminology, F^ are the a-cycles of TZ when j — 1, ... ,p — 2, while 
for j = p — 1, . . . , 2p — 4 they form the b-cycles. 
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Next, we introduce a notation for some special functions and meromorphic 
differentials on TZ. Function 



ivit) 1 



y A{t) t - z 

can be regarded as a meromorphic on TZ in both variables; 

duJk{t) ^ e{vk,t)dt, fc = l,...,p-2, 
is a basis of holomorphic differentials on TZ, and correspondingly, 



Uk(Z 



dujk 



form a basis of integrals of the first kind (these are multivalued and analytic func- 
tions on TZ having a constant increment along any cycle). 
Additionally, 



(3.16) 



IA{z) /-^ 



i{z,t)dt 



Viz) J a, 

can be also considered as an analytic function on TZ in both variables (multivalued 
in C,). Lemma 3.7 and Lemma 3.9 render that with an appropriate choice of z^^n = 



TT ^{zk,n), that means, either jz^^ or z)^'^, for 2: e C \ F, 

2 



.(2) 



p-2 p-2 

(3.17) NH^{z) = NG{z, (x.) + ^ 4,n Uk{z) + ^ ^ ^(2:s,,„, z) + 0{l/n), 



k=l 



fe=l 



where dfe.„ := N(vk.n — Vk) = NAvk.n, and as usual, we identify C \ F with the first 
sheet of TZ. 

Let us work out the system of equations on the unknown parameters. Given a 
closed contour (cycle) 7 on TZ, we denote by 



e(^;7) :=A^e{z,C) 



Ce7 



lA{z 



'V{t) dt 
A{t) t~ z 



, zeTZ\j, 



the period of 9{z, •) along 7. A direct verification shows that 0(-;7) can be ana- 
lytically continued on TZ as an integral of the first kind, so that for suitably chosen 

Cfc(7) e C, 



(3.18) 
By ( |3l2t . 



p-2 

0(2:; 7) ^^Ck{-/)Uk{z). 
k = l 



(3.19) N j) H'Jt)dt ^T{rj,f){l + 0{l/n)) mod 27ri, j = 1, 2, . . . , 2p - 4, 



and in view of (3.17), equation (3.19) may be written as 
(3.20) 

p-2 ^ p-2 

N 



i G'it,(x)dt = T{T,,f)-Y,dk,n i dL0k'~]^Y.^^^k.,^,Tj)+O{l/n) mod 2 

•'^1 fc=l •'^i k=l 
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with j ~ 1, 2, . . . , 2p — 4. This is a system of 2p — 4 equations on 2p — 4 unknowns 
di^n, • ■ • I dp_2^ny ^i.ni ■ • ■ : Zp-2,m that may be equivalently written in any basis Tj. 
From the general theory of Riemann surfaces it follows that matrix 

dojk I 

/ 3,k=l 



is invertible. Then, first p — 2 equations in (3.20) may be explicitly solved for dj, „. 
Substitution of those dk.nS in the remaining equations and the use of (3.18 1 reduces 
the situation to the standard Jacobi inversion problem, which as it is well known, 
is uniquely solvable for any non-special divisoij^ Hence, system (3.20) is uniquely 
solvable for any right hand side. 



Remark 3.11. Under Assumption 3.3 Avk.n — C>{l/n) and all the remainders in 
(3.20) are 0{l/n), which is the accuracy for determining dk^n by these equations. 



Now we can simplify the asymptotic formula (|3.4|) from Theorem 3.4 Since 
we get 



Ruiz) = a 



) (^)e-^^.W (l + e)(l/n)), 



and Hn can be replaced by Hn, the leading term in its asymptotic formula (3.17): 

^ P-2 ^ p-2 

(3.21) -Hniz) = G{z, oo) + -J2 dk,n Ukiz) + i^Y. ^(^fe^"' 

fe=l 



k=l 



Finally, using the analytic continuation of i?„ along a progressive path around 
an a € A (if we take a = ai, then 77„ just changes sign during this analytic 
continuation) and solving the system 

= Rn 
— Rn,l 

for Qn we obtain Q„ = (i?„.i — Rn)/ifi ^ /), from where the exterior asymptotics 
has the form 




Qn{z) = Cn,l 



I I'll 



-(z)e^^"(^) (l + e)(z)) 



/i/2(AV^)i/4' 

We summarize our findings in the following theorem: 

Theorem 3.12. For a normal index n N let N ^ ^ n(n + 1) and Hn he as 

and Zk.n determined by equations (3.20). 



defined in (|3.21), with the coefficients dk 
ippropriate 

Rn{z) — Cn,l 



Then, with an appropriate normalization 



and 
(3.22) 



i(z) = Cn.2 



{Avy/ 

1/2 



l_(^)e-^«.W (l + 0(l/n)), 
(z)e^«"W (l + 0(l/n)). 



/l/2(Ay)l/4 



In the situation when the divisor is special, we have deg hn < p — 2 and for such an n the 
normality is lost. 
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for z on a compact subsets of C\T. 

4. Case of p = 3: Chebotarev compact and the Riemann surface 

In the rest of the paper we concentrate on the particular case studied in 
when p = 3, ai, a2 and 03 are 3 non-colhnear points on the complex plane C, and T 
is the Chebotarev compact, i.e. the set of minimal capacity containing these points. 



Recall (see Theorem 2.1 ) that there exists a point v in the convex hull of A, called 



the center of the Chebotarev compact, such that with 

'Viz)^'/' 



A{z) = {z — ai){z — a2){z — 03), V{z)~z — v, and T{z) = 



A{z) 



where the branch of T in C \ F is specified by liuiz^ao zT{z) — 1, it is determined 
uniquely by the set of equations 



Furthermore, 
with 



Re / T{t) dt^Re / T{t) dt = 0. 



F = FiUF2UF3, 



i^zeC : Re J T{t)dt = o| 



the arc of F joining Oj with w, j ~ 1,2,3. We introduce also the orthogonal 
trajectories 



F-L :^ |z e C : Imy T(i) 0| , 
which consist of 3 unbounded rays emanating from v, as well as 

7^ := |z e C : Im^ T(i)dt = o| , j = l,2,3; 

each 7^ is an unbounded ray emanating from Oj , see Fig. [2| 

03 ^,--'73^ 



Figure 2. F and F-^. 
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Contour 7j'-urU7^ splits C\r into two simply connected domains. We denote 
by the domain containing on its boundary, and D_ the complementary one. 

On the three subarcs of T we fix the orientation "from aj to v" , while on the 
arcs of r-'- we choose the orientation "from v to infinity". This induces the left 
("+") and right ("— ") sides and boundary values. 

The equilibrium measure A = Ar on F has the form 

dX{z) = — T_{z)dz 
ni 



(compare with (2.5)). We denote also 



m,^\{Tj) = - I T_{t)dt, j = 1,2,3, 



SO that mi + m2 + m3 = 1. 
Define in C \ T 

(4.1) $(z) = exp T{t)dt 
normalized by the condition 

lim <i>(z) = 1; 

observe that $ coincides up to a multiplicative constant with exp(G(-,oo)) intro- 
duced in ( 2.4 1 . It is a conformal mapping of the exterior C \ F onto the exterior of 
the unit circle, such that 

(4.2) $(z) = cz + 0(l), z^oo, 

with 1/c coinciding, again up to a factor of absolute value 1, with the logarithmic 
capacity of F. Direct calculation allows to establish the following lemma: 

Lemma 4.1. For z G F° F \ {u, 01,02,03} and with the orientation shown on 
Figure 

(4.3) $_(z)<&+(z)-«;„ zeF°:-F,\{^^,a,}, 
with 

(4.4) Ki = e2"(™-'-"2), K2 = e-2"'"^ K3 = e2'^*'"^ 
so that = 1 and K2K3 — Ki- 

As before, we consider the Riemann surface TZ defined by the equation w'^ = 
A{z)V{z). Now it is an elliptic Riemann surface that can be considered as a two- 
sheeted covering of C, 7^ = {2; = iz,w) £ C^}, with two sheets, 7^^^^ and Tl'^^\ 
cut along F and glued together in the standard way. The canonical projection 
TT : 7?. — C is given by n{z) = z ioi z = {z, w) S C^. As in Section |3j we denote 
2;(j' = n^^{z)nTZ^-^\ j — 1, 2, and we convene that sheet TZ^^^ over C\F is specified 
by the condition w/z^ — > 1 as 2; — > cxj^^^ G TZ^^h In this way, function w = {AVY^'^ 
is single- valued on TZ, with 

't{z)/V{z), for 2 = ^(1) e 7^(l^ 
-T{z)/V{z), for 2 = ^(2) g -^(2)^ 

Again, we identify the first sheet TZ'^^^ with the domain C \ F = 7r(7?.'^^'). We also 
,0) 



w{z) 
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Figure 3. Cycles on TZ. 



We define the canonical homology basis of cycles as in Figure [Sj the a-cycle 
encloses v and ai, while the b-cycle goes around v and 03. Both are oriented as 
indicated in the figure, so that at their unique intersection point on TZ their tangent 
vectors form a right pair. 

The normal form of the Riemann surface TZ is the polygon (rectangle) TZ with 
sides aboT^h^^ (see Figure |4]j^ 

0.3 




Figure 4. Polygon TZ. Shaded domain corresponds to the first 
sheet 7?.^^\ and the dotted line corresponds to 7r^^(7;j'- U 7^). 

We introduce also some notation, slightly different from that used in Section 
|3j related to differentials on TZ and their periods. For integer fc, denote 

(4.5) dv,{z) = 

w{z) 



^Onc of the authors of [2j kindly pointed out to us that a similar figure is contained in the 
cited paper. 
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Then di^o is, up to a constant multiple, the only holomorphic differential (abelian 
differential of the first kind) on TZ. 

For two points ri,r2 € TZ we denote by fi^i ,r2 the normalized differential of 
the third kind, such that it has only simple poles at ri, with residue +1, and at r2, 
with residue —1, and with a vanishing b-period. In particular, if z* = (z* , w*) €z TZ, 

1 f w{z) + w* 



dfl^, (1) = ^ z + S] dz 

2w[z) \ z — z* J 



b w{z) \ z- z* J / / V/fa w{z) 

With the orientation of Tj specified in Figure [2] we define 

(4.7) M^--=-7^- / ^TV, fceNU{0}, J = 1,2,3, 
so that 

i dvk = -47ri M^^, i dt/^ = -47ri Mg*^, fc e N U {0}. 
In particular, Af" 7^ 0, j ' = 1, 2, 3, and Im(r) < 0, with 

(4.8) r := ^. 

Direct computation using the Cauchy integral formula shows that 

Af + A/2" + A/3" - 0, A/1 + A/1 + A/1 = ^ , A/2 + A/| + A/| = | , 
where we use the notation 

(4.9) S " + °^+"^+"^ 

We reserve the notation di^Q for the normalized differential of the first kind, 
whose b-period is equal 27ri: 

'^"(^)--2if^ = -2]^'^^"- 
Observe that the a-period of i'q is 27riT. 

5. Riemann and Hilbert, or the non-linear steepest descent analysis 

Now we are ready to return to the Fade approximants of the function 

f{z)^{z-air^iz~a2nz-a,r^, 

with ai, a2, ck3 G M \ Z, such that ai + ^2 + 0^3 = 0. For the sake of simplicity of 
the analysis we assume additionally that aj > — 1, j = 1, 2, 3. As before, we specify 
the branch in C \ F by /(oo) = 1 and agree in denoting by /^/^ the branch of the 
square root in C \ F given by /^/^(oo) = 1. 

We need to introduce an additional piece of notation: for j — 1, 2, 3, let 

(5.1) Tj = e"*'^"^', tj = 2i sin(7rQ!j), s„j = tjn", 
with Kj defined in (4.4). Observe that T1T2T3 — 1, and 

(5.2) Tj + ltj^l + tj + tj + lTj\ = + 1 + tj + tj-lTj^^ =0, j = 1, 2, 3, 

where the subindices are taken mod 3. 
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Let us collapse the contour of integration in (2.3) onto F; as a consequence, 
function / induces on F the weight p, 

p(z) = /_(2)-/+(z) = ^/+(z) = (rr2_i)/+(^), zer,\{v,a,}, j = 1,2,3, 

so that the orthogonality condition ( |2.3| can be rewritten as 

(5.3) J t^Qn{t)p{t)dt^G, fc==0, l,...,n-l. 

By our assumption that aj > —1, the weight is integrable on F, and the regularity 
of / at infinity implies that 

lim p{z) + lim p{z) + lim p{z) — 0. 

zGTi zeT2 zeTi 

Standard arguments show that there is an integral formula for the residue _R„: 

Rn{z) = — / — 7 dt, z e C \ F. 

Ztti Jy t — z 

This allows us to formulate the Riemann-Hilbert problem for Qn and Let 
(73 denote the third Pauli matrix, 



^3 = 



1 

-1 



and for any scalar a we use the notation 



We seek the matrix-valued and analytic function Y = Y{-; n) : C \ F ^ C^^^, such 
that: 

(RH-Yl) It has continuous boundary values Y± on both sides of F°, and with the 
specified orientation of F, 

F+(z) = r_(z) (^J P^^'^y zeT°. 

(RH-Y2) Y{z) = {I + 0{l/z)) z'^^S as z ^ oo. 
(RH-Y3) As z ^ flj, z e C \ F, j = 1, 2, 3, 



Y{z)- 

(RH-Y4) As z ^ -y, z e C \ F, 



^ 1 l^JUh if".<0' 



1 


\z 




1 






o 












if aj > 0. 



Y{z)=o(] I^SI---! 
^ ' \l log \z — v\ 

From the fundamental work of Fokas, Its and Kitaev |14| it follows that 

Theorem 5.1. The matrix valued function Y{z) given by 

Quiz) R^{z) \ 

,-2Tri"fl_iQn-liz) -2TTijl_^Rn-liz)J 



Yiz) = 
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is the unique solution of (RH-Y1)-(RH-Y4-), where Qn is the monic polynomial 
of degree n satisfying (5.3) and 7„ is the leading coefficient of the corresponding 
orthonormal polynomial. 

This result is complemented with the non-linear steepest descent method of 
Deift and Zhou [6l |9l I10|, 1111, 112j : we need to perform a number of explicit and 
invertible transformations of (RH-Y1)-(RH-Y4) in order to reach a boundary value 
problem with jumps asymptotically close to the identity and a regular behavior 
at infinity. Two of the main ingredients of this analysis are the outer (global) 
parametrix and the local model at the Chebotarev center, that we explain next. 

5.1. Global parametrix. For n e N we need to find an analytic matrix- 
valued function = N : C\r^C2^^ such that 

(RH-Nl) It has continuous boundary values N± on both sides of r°, and with the 
orientation "from aj to u" of F, 

(5.4) iV,(.)^Ar_(.)(_J^^^^^. ^-), .eF°. 

(RH-N2) N{z) = 1 + 0{l/z)), as z^oo. 
(RH-N3) As z ^ Oj, z eC\ F, j = 1, 2, 3, 

Ar(z)=0(|z-a,|-V4). 
As z ^v, z eC\T, 

N{z) =0{\z-v\-^^^). 



Constants s„j- were defined in (5.1); hence, the dependence on n resides only in 
the boundary condition (5.4 1. 

On the compact subsets of C \ F this problem is asymptotically close to the 
boundary value problem for the following matrix, 
(5.5) T{z) :=c"''^r(z)$-"'^^(z)r^/2(^-)^ 



with c defined in (4.2). Hence, we can expect that away from the Chebotarev 
compact F the solution N of (RH-N1)-(RH-N3) models the behavior of T for n 
large enough. 

We build TV in the following form, 

(5.6) N{z) := Fiooy-'N{z)F{z)-''\ 

using two "ingredients" described in detail below: a scalar function F, which plays 
the role of a Szego function with piece- wise constant boundary values, and a matrix- 
valued function TV = (Nij), which will be defined in terms of abelian integrals on 
TZ. Both F and TV depend on n, but in this section we omit this dependence from 
the notation, keeping it in mind. 

With the notation (4.51 consider the equation 

di'o = - ( 1 + log f — ^ ^ ^ d.i'o- ^ log f — ^ / di^o, 



27ri \Sn.lJ J In 27ri VSn,3 



or equivalently, 

* ■ f 

rli^r. = — 'y.'Tvi.T I 1 -I- 

2711 \SnAjJ \Sn,3 



(5.7) / d,.; = -27r^r I 1 + log f '^') ) - log ^'"-^ 
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where i/q is the normahzed differential of the first kind (4.101, and the path of 
integration lies entirely in the rectangle TZ. Among all possible choices of the 
branch of the logarithm, there is at most one value of log(s„^2/sn,i) and at most 
one value of log(s„^2/sn.3) such that this equation has a solution in TZ; this solution 
Zn = {zn, Wn) is obviously unique. 

Remark 5.2. If Zn falls on one of the cycles Tj we consider it slightly deformed so 
that the same argument applies. A truly special situation occurs when eventually 
Zn = oo'^^ or Zn = oo^^^. The first case happens when 

(5.8) |S,„,f£-U|«l„gf5-Uo „,„<iZ, 



27ri \Sn,iJ 27ri \Sn,3 



The consequences of this degeneration are discussed below, see Remark 5.8 



With this choice of the branch of the value of log(s„^2/sn,i) we define two 
parameters, /3i and /32, as follows: 



(5.9) /3i := log Sn,i = log(ti) + 27rm(m3 - 1712). 

(5.10) /32:=7ri + log 



Sn,2 
Sn,l 



Obviously, (5.9) defines /3i up to an additive constant which is an integer multiple 
of 27ri. 

With these two complex constants fixed, we build a complex-valued function 
-F on C \ r, holomorphic, uniformly bounded and non-vanishing in C \ F, and such 
that 

z e Ti \ {v,ai}, 

(5.11) F+{z)F^{z) = Uf^+f- ^ -sn,2, zer2\{v,a2}. 




Constant /^s is not arbitrary: 

(5.12) /33-(1 + t)/32 



with T from (4.8|. We take F of the form F{z) = exp(A(z)), and give two equivalent 
expressions for A. 

First, A can be built in terms of the holomorphic differential 1^0 on M: 

(5.13) A(z) = | + s r"dz.o-| + s r"4Y' ^'^^e^ew, z€C\r, 

2 7ai 2 J^^ w{t) 

where 

(5.14) S = /?2 [(Mf - rM|) - S {Ml - tM^)] , 



with Mf introduced in (|4.7| and S in (|4.9[). The path of integration in (|5.13p 



lies entirely in TZ^-^\ except for its initial point. Observe that A in (5.13) is a 
holomorphic function in C \ F. 

Alternatively, define the functions 
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where we integrate in the direction "from Uj to f " , and let 



(5.16) 



A(z) 



f32{h{z)-Ti3{z)-i/2), zec\r. 



Lemma 5.3. With A given either by (5.13) -(5.14) or by (5.16), function F{z) = 
Fn(z) = exp(A(z)) is holomorphic, uniformly bounded and non-vanishing in C\T, 
with 



(5.17) 



F(oo 



)=exp(^^-/32(Mi-TA/i-l/2)). 



Moreover, F has continuous boundary values at T° that satisfy (5.11)-(5.12|. 

In consequence, formulas (5.13)-(5.14) and (5.16) define the same function in 

c\r. 

Remark 5.4. Recall that /32 was defined uniquely as a function of n, but /3i is 
determined mod (27ri). From (5.13)-(5.14) or (5.161 it follows that for each n, 
function F is determined uniquely up to a change of sign. 

Now we define in C \ F the analytic matrix valued function TV entry-wise in 
terms of meromorphic differentials on TZ as follows (see also |19j l. The meromorphic 
differential 



drf(z)--'^-^^dz^'d.* 



1 



^ 1 



1 



a, Mlwiz) 



dz 



has only simple poles on TZ: at the zeros of AV with residues —1/2, and at oo(^\ 
oo^^-*, both with residues +1; additionally, its b-period is zero. 

With Zn — {zn,Wn) solving (5.7) we consider the meromorphic differential 

(5.18) 

or more explicitly, 
(5.19) = 



(AVYiz) 



Si 



dz. 



4{AV){z) 2(z-z„) 2w{z){z-Zn) 2w{z) w{z) ^ 

where Si is uniquely determined by the condition that the b-period of rjz^ is zero. 
Obviously, it has the only poles, all simple, at the zeros of AV with residues —1/2, 
and at z„ and oo'^^\ both with residues -fl. 

Lemma 5.5. With the conditions above, 

(5.20) / d7fe„ = /3i + /33 - log(s„,3) mod (2^z), 



with jSi and 13^ given in (5.9) and (5.12), respectively. 
Proof. Direct calculation shows that 

drj* — ni {1 ~ t) , 



and from the Riemann's identities it follows that 

1 



(5.21) 



dvX 



It remains to use (5.7 1 and the definition of /?i and ^2 above. 



□ 
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Remark 5.6. The uniqueness of Zn satisfying (5.20) can be easily established 



for any other r £ TZ, v^^ — t'^ is a meromorphic differential in 72. whose only 
poles (both simple) are at 2:„ (with residue 1) and r (with residue —1), and with 
periods multiple of 27ri, so that exp(J diy^,^ — v,.)) has a single pole at r, which is 
impossible. 



Let 
(5.22) 



ui{z) — exp 



U2(z) = exp 



with z^-''' = TT ^{z) n TZ^j\ For ui, the path of integration lies entirely in TZ^^\ 
while for U2 it goes from into crossing r2 once, see Figure [H] 




Figure 5. Paths of intcrration for functions Uj defined in (5.22) 



Lemma 5.7. Functions Uj are holomorphic in C \ {j^ U F U 72 ) (see Figure^, 

and have a continuous boundary values on •y^ U F" U 7^ such that 

(5.23) 



(wi)_(2), zej^u-f^, 



-{u2)-iz), ze7^U7f 



K)+W = <|("li2' {u2hiz) = \{mUz), zenuF^, 



Sn,3 



-(W2)_(z), zeF°; 



S n,3 
l+ft 



(wi)_(z), zeri. 



Moreover, as z aj, z G C \ F, j = 1,2,3, Uk{z) = 0{\z — aj\ ^^'^), while as 
z — >■ u, z € C \ F, Uk{z) — 0{\z — wj^"'"/*). Additionally, 

ui(z) = 1 + 0^ , U2(z) = O0^, z^oo. 

Finally, if Zn e TZ'^^^ then ui has a simple zero at z = Zn and M2(z„) ^ 0; otherwise, 
U2 has a simple zero at z = z^ and Mi(z„) ^ 0. 

With these two functions we define in C \ F 

W2(z), if z G Z)_|_, 



(5.24) 



Nii(z) = wi(z), Ni2(z) 



—U2{z), if z G _D_. 
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Recall that the simply connected domains D± are limited by 7 j*" U F U 7^ , and 
is the one containing 03 on its boundary. 

Remark 5.8. As we have noticed before, it may happen that either Zn — oo^^^^ or 
Zn = oc^'^K In the first case, when condition (5.81 holds, we just have rjz^ = V* : 
so that TVii has a zero at infinity. This, as it follows from the asymptotic formulas 
below, will mean that the index n is not normal, see the expression of x in (6.11. 
In the second case, r]^^ = 7700(2) has a simple pole at cxd^^-' with residue +2, which 
creates a double zero of N12 at infinity. 

Furthermore, consider a family of functions g on 7?. of the form 

'w{z) + Wn 



a,b £ 



Each such a function has a simple pole at Zn = (z„,w„) and at (x>^'^\ There is a 
unique combination of constants a, b, such that additionally 

g(oo(i)) = and lim q{z)Ni2{z) = 1. 
Let q^^\z) — q(zS-^^) be the values of q on the j-th sheet. Then set 



This defines completely the matrix- valued function TV = (Nij). Finally, we assem- 
bly TV as in (5.61 using this matrix TV and function F given by ( 5.13 1-( 5.14) or 
(5.161 with parameters (5.9)~(5.10). Direct verification shows that the following 
statement holds true: 

Proposition 5.9. Matrix TV constructed above solves the RH problem (RH-Nl)- 
(RH-N3). 



5.2. Local parametrix. Matrix T defined in ( 5.5 1 has jumps that are asymp- 
totically close to the identity matrix for n large enough, as long as we stay away 
from F. However, this behavior fails in a neighborhood of A and the Chebotarev 
center v, where we need to perform a separate analysis in order to find an appro- 
priate model. Here we describe only the construction of the local parametrix P at 
z = V (around the branch points Oj matrix P is built in the way described in detail 
inm)- 

We take a small 6 > and define Ds := {z G C : \z — v\ < S}, Bs {z E C : 
\z — v\ =6}, assuming that Ds (1 A = 0, see Figure [6] 
The local parametrix at z — v has the form 



(5.25) 



P(z) := E{z)'^ (n^/^(piz)^ b"'' B 'S?-'^"^ (z), z e Ds\{TUT-^). 



Here $ is the conformal mapping defined in (4.1), and the rest of the ingredients 
are: 

• Function ip, defined piece- wise in each sector: for z G Ds \ (F U F-*"), 



(5.26) 



ip{z) := 




T{t)dt 



2/3 



T{t)dt] 



2/3 



if z e ® u ' 

otherwise, 
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where we take the main branch of the power function. Then (ys is a con- 
formal mapping of Dg onto a neighborhood of the origin, ip{v) = 0, and 
Tj^ is mapped onto the positive semi axis. 
• Constants bi and 62, determined up to a change of sign by 



&16 



and 



• Matrices 



1^2 



b := 



lAi, 61/&2 



-^2/^3 



61, ifze®u®u®, 

62, ifze@u@u®. 



B ■= < 



-i, 

1^2 

i3' 
-IA3 



ifze®u@u@U' 
if z e ®, 

if z e ®, 



and 



M2 = M2{z,n) 



('^?/&2)'^^ 
' 62*2' 
-l/(&2i2) ^ 



if z e ® u ®, 
if z e @ u @, 

if z e z e ® u ®. 



• Matrix valued functions in Ds \ Ti, 



Mi(z) := 



(73/4 



(5.27) 



g-7ri/3 g-7ri/3 

where we take the main branch of the root, and 

E{z) = £;(z, n) 7V(z)M2(z)Mi(z)n'"^/^ z e \ T, 



with N constructed in subsection 5.1 see (5.6). E{z) extends in fact as 
a holomorphic function to the whole Dg. 
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• The Airy parametrix, '5', defined as 
'Az{0 Ai{uj\) ■ 



*(C) 
*(C) 
*(C) 
*(C) 



-(73/4 



-0-3/4 



Ai'iC) oj^Ai'{Lo^C) 

Az{0 Ai{uj\) 

Ai'(C) uj^Ai'iuj^C) 

Ai{0 -a;2yli(wC) 

Ai'iO -M'{ujO 

Ai{C,) ^uj'^Ai{ujQ 
where uo — exp(27ri/3), see e.g. [71 



1 

-1 1 



0-3/4 



1 
1 1 



Theorem 5.10. Matrix-valued function P given by (5.25) solves the following 
boundary-value problem in Dg: 

(RH-Pl) It has continuous boundary values P± on both sides of all curves, and with 
the specified orientation, 



P+{z)^p^{z)jp{z), ze (r°ur^)ni?5, 



where 



(5.28) 



(5.29) 



Jp 




(RH-P2) P{z) = + 0{l/n)) N(z), for z G Bg. 
(RH-P3) As z V, z eC\{ru T^), P{z) = 0(1). 

5.3. Asymptotic analysis. In the final transformation of the original prob- 
lem (RH-Y1)-(RH-Y4) we define a matrix valued function R in the form T{z)A~^{z), 
where, roughly speaking, A = N away from F and A = P in a neighborhood of A 
and V. The explicit formula for A in Dg is given above, while it is built in terms of 
the Bessel functions in a neighborhood of the branch points Oj, see |18| for details. 
The inverses of all these matrices exist, since the determinants of these matrices 
are equal to 1. 

The construction of TV and P is such that 



R+{z) = R^{z) [I + O 



uniformly on a finite set of contours in ' 
we conclude that 



Following the already standard reasoning 



(5.30) 



R{z) 



O 



00, 



uniformly in C. The relation (5.301 is the main term in the asymptotics for R and 
it is enough to give the leading term in the asymptotics for Y . 

For instance, taking into account (5.5 1 and (5.30) we see that locally uniformly 
in C \ F we have 



Y(z) 



-n^3 + J^l^^ Ar(z)/(z)-'"3/2$"'^3(^)^ 
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where c was defined in (4.2 1. In particular, working out the expressions for Yn and 
Yi2 we get: 

Theorem 5.11. Locally uniformly in C \ 
(5.31) 

F(oo) 



Quiz) = 

and 

Rn{z) 



exp 



l + O 



,1/2 



icHz)) 



j,F{^)F(z) 



± exp 



l + O I - 



O 



n 



The sign and the paths of integration are selected in accordance with the defi- 
nition of Nij in (|5.24). 



We see in particular, that the spurious zero of Qn is asymptotically close to the 
unique zero of iVn, which appears only when Zn is on the first sheet. Otherwise, 
it gives us an extra interpolation condition (zero of 1^2, close to the zero of Ni2)- 

The Riemann-Hilbert analysis yields asymptotic formulas not only away from 
r but in the rest of the regions. For instance, close to F but still away from the 
branch points aj and the Chebotarev center v the asymptotic expression for Qn is a 
combination of two competing terms, which gives rise to zeros of Qn- For instance. 



by (5.5) and (5.30), for z in ® of the domain D, 



Y{z) 



= c-'-^ [l + O ) N{z) J) /(z)--^'/2$"-3 (^)_ 

In particular, 

c" Qn{z)f{z)''^ = (^Arn(z)<f>"(z) + Ari2(z)^$-"(z)) {l + O 

+ ('Ar2i(z)(f>"(z)+iV22(z)^$-"(z)') O Q 



For z e Fi we can rewrite it as 

'$+(z)V / N^i{z)+ , Arii(z)_ 



Qn{z) = 



,1/2 



1/2 



l + O 



Finally, in order to analyze the behavior at the Chebotarev center v, for in- 
stance, when z £ Dgn® (see Figure [6| , it is sufficient to obtain the expression for 
Y from 



J) /(2)-'>'-«""(i), 



where all the ingredients in the right hand side were given above. Since the formulas 
for Qn and i?„ obtained this way are not easily simplified, we omit their explicit 
calculation here for the sake of brevity. 
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6. Wrapping up, or matching the asymptotic formulas and the Nutall's 

conjecture 

In |27j Nuttall conjectured the form of the leading term of asymptotics for Qn 
and Rn away from T in terms of a solution of a scalar boundary value problem. 
We show next that our results match the Nuttall's conjecture. 

Let us denote 

7Vn(z) 



(6.1) xi-A 



so that X has a pole of order n at infinity, and 5H has there a zero of order n + 1 
;ical situation of z„ ~ < 

Qn{z)=x{z) (l+O 



(unless the pathological situation of z„ ~ oo occurs). By (5.311, 

1 



z e c\r. 



But for z e r°. 



(/(z)x(z))± = /(z)f (^^M±y N^,^{zy, 



using that p{z) ~ tjf+{z)/Tj ~ tjTjf^{z) on r°, (5.4) and Lemma 4.1 we get 



'^{z)x+iz) = ~ w+{z) 



f(z)^^^ 

^-^-^ {c<i>{z)^)-" N,2-{z) - {wD\)-{z), 



a{z)x-{z) =w+{z)T,f{zfJ' {cMz)+r" N,2+{z) = {wm)+{z), 

where cr(z) := p{z)w^{z) on r°. These two equations match the boundary value 
conditions in [28', formula (5.6)] (after replacing X2 = X and H = w^i) that define 
uniquely the leading asymptotic terms for Qn and i?„, according to the conjecture 
of Nuttall. 

Let us finally compare the asymptotic formulas obtained in Sections |3] and [5] 
and given by Theorems 3.12| and 5.11 We introduce here the notation 

dt 



u{z) 



On one hand, observe that in the case p = 3, with the function 0{z, C,) defined in 

dt 



(3.16) we have 



w„ 



V{zn 



-u{z) + w„ 



ai W{t){t~ ZnY 



SO that 



NUniz) ^NG{z, oo) + d„ u{z) + - e{zn, z) 



=nlog $(z) 



dt 



where 



tdt 



2 Ja, w{t){t-zn) 2 w{t) 



^2 = dn - 



+ 52U{Z) + O - 



2 ■ 21/(z„)- 

On the other hand, up to a multiplicative constant. 



F{z) ^ exp / di-jz 



exp 
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Recalling the expression for 77^^ in (5.19), we get 



F(z)-iexp (^J 




with an appropriate selection of the constant 5^. Comparing thus expressions 



(3.21 H(3.22 1, obtained by the WKB analysis, with (5.31) we see that they co- 
incide, up to the right determination of the constants 82 and 5z above. But these 
constants are uniquely determined by the condition that the right hand side in 



(3.22 ) and in (5.31 ) must be single- valued in C \ F. 
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